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We study a Randall-Sundrum model modified by a Gauss-Bonnet interaction term. We consider, 
in particular, a Friedmann-Robertson- Walker metric on the brane and analyse the resulting cosmo- 
logical scenario. It is shown that the usual Friedmann equations are recovered on the brane. The 
equation of state relating the enery density and the pressure is uniquely determined by the matching 
conditions. A cosmological solution with negative pressure is found. 

October 2001 
I. INTRODUCTION 

The possibility that our universe is a four dimensional brane embedded in a higher dimensional spacetime has been 
extensively discussed recently. The most popular model in this context is the one proposed by Randall and Sundrum 
. This scenario is based on the metric 

ds^ = {y) riijdx'dx' + dy'^ , (1.1) 



T-H , where r]ij is a flat Minkowski four dimensional metric on the brane and A'^{y) is the warp factor depending only on \y\. 

^ ■ Perturbations of this metric reproduces the expected 1/r Newtonian potential on the brane (the observed universe) . 

' This is due to the fact that the zero modes of the perturbation propagate on the brane only (they tend rapidely to 

^ ' zero in the fifth dimension). The other modes (the massive Kaluza-Klein modes) give merely a correction in 1/r^ to 

' this potential. 

I t One of the first development of this model was the generalisation of the Randall-Sundrum ansatz to include a wider 

y—i ' class of metrics P-[l5||. Different geometries were treated by considering solutions to the Randall-Sundrum model 
with metrics which, up to a non-constant conformal factor, can be written as 

; ds^ = (y) g^j {x, y) dx'dx^ + dy^ . (1.2) 

The requirement that the zero modes of the perturbations around these metrics are localised on the brane imposes 
(— I I further constraints on this class of geometries 

• • • Since the Randall-Sundrum model is a string inspired picture []l6| , one would like to understand the implications of 
. !^ higher curvature terms in such brane world universe. These terms naturally arise in the string effective action beyond 
. the first order in the string tension a' . The inclusion of these terms are also of relevance to cosmology and inflation. 
^ ' It turns out that generic higher curvature terms lead in general to a delocalisation of gravity from the brane ||l^,|l8| . 
. 5^ , A combination of these terms in a Gauss-Bonnet form yields, however, the desired Randall-Sundrum behaviour of 
the zero modes of the perturbations. We sho uld mention that the analyses of refs. |l^,|l^ is carried out with a brane 



possessing a flat metric of the form given in (1.1). In this context, various other issues were also treated using higher 
curvature terms ||l^~|27||. 

In this letter, we deal with a Ran dall- Sundrum model complemented by a Gauss-Bonnet density where the five 



dimensional metric is of the form in (1.2). We start by considering a metric on the brane with spherical symmetry. 
It is shown that the only possible solution in this case is a de Sitter or anti de Sitter spacetime for the brane. This 
is in contrast to the case without a Gauss-Bonnet term where black hole geometries are allowed Our bulk metric 
involves a warp factor that presents an oscillatory regime among other possibilities. 

A second study consists in taking a Friedmann-Robertson- Walker metric on the brane. We recover the equations of 
ordinary cosmology on the brane. This is to be compared to previous brane world cosmology models [p8[-^8| where, 
among other things, the square of the Hubble parameter is found to be proportional to the square of the energy density. 



*Permanent address: Laboratoire de Physique Theorique, Centre Universitaire Mustapha Stambouli, 29000 Mascara, Algeria, 
t E-mail: boucif@celfi.phys.univ-tours.fr 
■'•E-mail: nouri@celfi.phys.univ-tours.fr 



1 



Here, the matching conditions are so restrictive that they determine the equation of state relating the energy density 
to the pressure. Various inflationary solutions with a cosmological constant on the brane are determined. Another 
solution with a time dependent energy density and pressure is also found. However, the pressure for this matter is 
negative and cannot describe an ordinary dust. On the other hand, a scalar field is found whose energy-momentum 
tensor could describe this behaviour. 

II. THE MODEL 

We consider a five dimensional spacetime with coordinates = t,x^,x'^,x^,x'^ = y) where {t,x^,x'^,x^) denotes 
the usual four-dimensional spacetime and x^ = y is the coordinate of the fifth dimension, which is an orbifold jZi 
where the Zi action identifies y and —y. The five dimensional indices are denoted by M, N, . . . — 0, . . . , 4 and the 
four dimensional brane world indices are i, j, . . . = 0, . . . , 3. We will neglect the matter interaction and consider the 
five dimensional gravitational action 

S ^ J d^xV^{aR + A + (3Lgb) (2.1) 

where a and (3 are two coupling constants and A is the five dimensional cosmological constant. The Gauss-Bonnet 
Lagrangian density is 

Lqb = RmnpqR^^^^^ — ^RmnR^^'^ + R^ ■ (2.2) 
The equations of motion corresponding to our action ar^ 

^MN = a ^Rmn — -j^gMNR^ — -^AgMN 

+ 2p^Rj^^'^^ Rnpqs + 2R^'^ Rmpqn — '^RmpRJ' + RRmn — ^QmnLq^^ —0 . (2.3) 

As it is well-known, these equations possess the following solutions p9| , |50t| 

Rmnpq = {GmpGnq — GmqGnp) , (2.4) 
where a takes two possible values as given by 



.= i^(^a±f^ J/?Aj . (2.5) 

This last equation can be inverted to get an expression for A 

A 12cr (a - 2/3cr) . (2.6) 

This value of A will be needed in the rest of the paper. 

As a start, we would like to analyse the above equations of motion in the context of the brane world scenario. We 
consider, for this purpose, a spherically symmetric line element as given by 

ds^ = A^{y) (^~N{r)dt^ + M{r)dr^ + r^(de^ + sin^ Odcj)'^)^ + dy^ . (2.7) 

and examin the equations of motion 

n'j^ = £ll +T^' = Q , (2.8) 
where the non-zero components of the energy-momentum tensor T^^ are 



^Our conventions are such that R^'mpq = dpT^g + Tp^r^g ~ {P ^ Q) and Rmn = -R'^j 



MQN- 
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T;^S}XS{y) (2.9) 

with A denoting the cosmological constant on the brane. We obtain four different equations corresponding to the 
components Hq, Hi, = H.^ and 7i|. Substracting Hi from Hq leads to 

M(r) = ^. (2.10) 

Substituting for M (r) and substracting H2 from Hq yields 

N{r) = l + fir^ + - , (2.11) 
r 

where fj, and v are two constants of integration. 

After substituting for N (r) , the equation corresponding to Til — can be cast in the form 

12 (^ + A'2) [aA^ - 2/3 + A'^)] - AA^ = 1^ . (2.12) 

Our notation is explained in the footnote below^ . It is clear here that one must have v — 0. This condition, however, 
is not needed if /? = (see for example Q). With h' — 0, the last equation takes then the simple form 

H + A'^= (jA^ , (2.13) 

where a is as previously defined. The solution to our last differential equation is 

^(2/) = ^exp(±V^|y|) +7exp(TV^|y|) . (2.14) 

Here 7 is an integration constant. 

We are left with one equation to solve, namely Hq 

QaA {fi + A'^) + 6 [aA^ - 4/3 (^ + A'^)] ^ AA^ + 2XA^5 (y) = . (2.15) 

ay 

This equation involves second derivatives of A (y) which generate delta functions as explained in the footnote. Using 
equation ( ^.13| ) in Hq and matching the delta functions yields the fine tuning conditions 

A = -6 (a -4/3.7)^ . (2.16) 



Substituting then for A (y) in ( p. 15 ) , fixes the five dimensional cosmological constant A to its original value A = 



12cr (a — 2(3a). Notice that the warp factor A{y) can present various behaviours depending on the value of a. In 
particular if a is negative then an oscillatory regime is obtained. 

III. COSMOLOGICAL SOLUTIONS 

The metric for this study is taken to have the form 



dr 



2 



1 — kr'^ 



+ r^(dr +sin^ Odct)^) 



dy' . (3.1) 



The non-vanishing components of the energy-momentum tensor Tj^ are now given by 



Since A is a function of \y\ we have ^ — A'^^ where A' denotes the derivative of A with respect to its argument \y\ and 
= 20 (j/) — 1 where (j/) is the Heaviside function. Notice that (41^) = 1 B.nd we have (^) ~ A'^. On the other hand 
= A" + 2A'S (y) where A" denotes the second derivative of A with respect to its argument |y|. 
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(3.2) 



The equations of motion we would like to solve are still Hj^ = £j^^ + Tj!f = 0. There are three different types of 
equations Hq, Hi —H2 = and Hf. The first of these, Hq, is 

6 [a^ {aA^ - W^) + 4/3 (fc + a^)] ^ - 6aA {k + a^) + Qaa^AA'^ - AA^a^ - 2pa^A^S (y) = , (3.3) 

ay 

where a is the derivative of a with respect to t. Matching the delta functions in this last equation yields 

6A'{0) 



a^A{0)f 

The second equation, 7i}, is given by the expression 



(a [^(0)]' - 4/3 [^'(0)]') + 4/? (fc + a^) 



(3.4) 



d^A 



2 [Sa^ {aA^ - 4,(3A'^) + 8/3aa + 4/3 (A: + d^)] _ 2aA {k + a^) + Gaa'^AA''^ ~ Aa^A^ - AaadA 
+2pa'^A^5 (y) = . 



(3.5) 



Again, matching the delta functions coming from and those coming from the energy-momentum tensor gives 



3a2 (a [A{Q)f - 4/3 [A'{Q)f) + 8/3aa + 4/3 (/c + d^) 



_ 2A'{fd) 
^~ a^A{Q)f 

Using the expression of p we deduce the following expression for the Hubble parameter 



[Am 

24/3/1' (O)"" a 



(3.6) 



(3.7) 



This is the first Friedmann relation of ordinary cosmology. Similarly, combining the expression of p with that corre- 
sponding to p, yields the second Friedmann equation 



(3.8) 



Differentiating and using the expression of a results in the usual conservation equations 

pa -I- 3 (p + p) d = . 



(3.9) 



Therefore, our gravitational theory with the Gauss-Bonnet term leads to ordinary Friedmann equations. 

Once the matching is carried out, we can deal with the two equations 7i} and Tip away from y = 0. Substracting 
"HX from T^g gives then 



{k + (? - aa) {aA - Af3A") = . 
If the first factor (fc + d^ — ad) vanishes then an interesting solution is given by 

k 



a it) 



4kt2 



exp (±rt) -|- Kcxp (=pri) , 



where t and k are two integration constants. 
Upon substituting for a{t) in Til, we obtain 



12 {A'^ ~ T^) [aA^ - 2/3 {A'^ - t^)] - AA* = 



(3.10) 



(3.11) 



(3.12) 



This equation is exactly the one found in (2.12) with = Q and where p is replaced by — r^. The solution, A{y), to 
this equation is therefore as given in (2.14) upon replacing p by — r^. 
Putting the expression of a{t) in Hq leads to the differential equation 
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6aA - t2) + 6 [aA^ - 4/3 {A'^ - r^)] A" - AA^ = 



(3.13) 



Again, this equation is that found in (2.15) away from the position of the brane and where is replaced by — r^. There- 
fore, replacing for A{y) in this last differential equation fixes the bulk cosmological constant to be A = 12(7 {a — 2j3a). 
We should mention that this solution leads to the following relation between the energy density p and the pressure 

P 



P 



X 



where A is as given in ( 2.16 ). The Hubble parameter for this solution is given by 



with a (t) as given in (3.11). Similarly, we find that 



a 
a 



(3.14) 



(3.15) 



(3.16) 



We distinguish, therefore, two cases. The first corresponds to > and leads to an inflationary regime whenever 
one of the two exponentials in a [t] dominates. The second situation arises when < 0. In this case the scale factor 
is given by 



a{t) = s cos 



-t + if 



(3.17) 



where e and ip are two real integration constants. The corresponding A{y) is given by equation (2.14) where /i is 
replaced by —k/e^. Of course, this solution is valid only when k is negative and describes a repeatedly collapsing 
universe. The Hubble parameter in this case is = — fc/a^ + k/e"^. 



The other solution to (A: 



aa) = is given by 



a (t) = ±V^t + S 



(3.18) 



with S an integration constants. This solution is physical only for k negative. Substituting this expression of a {t) in 
leads to the differential equation 

A'^ = aA^ , (3.19) 



(3.20) 



where a is as defined in ( |2.5D . Therefore A (y) is given by 

A (y) ^ ^ exp {^/a\y\) or A (y) = uj exp ^/a\y\) 



for two integration constants ip ^^^d The solution in ( 3.18| ) and ( 3.20| ) automatically satisfies Hq. Furthermore, we 

have p — —p — A and the Hubble parameter is — —k/a'^. 

Let us now return to the second possibility as allowed by equation (3.10), namely when [aA — Aj3A") — 0. This 
case is of course present only when f3 is different from zero and we have the solution 



■f fa 
A{y) =Cexp ( 2A/^lyl 



6 exp 



1 



y\ 



2V /?' 

with 1^ and 9 are two integration constants. Substituting this solution in leads to 



(3.21) 



(Sa^ - 2/3A) 



Cexp 



1 r^. 



\y\ 



exp 



1 fa, 



\y\ 



2q2 



6^0 



-48a [P^ (fc + d^) + af3^0a^] - 6^0a [8a/3 {k + a^) + iOa^ + 2/3A)] = 
It is clear, from the separation of variables, that one must have 

3a2 



A: 



2(3 



(3.22) 



(3.23) 
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and the above differential equation reduces then to 

[P {k + a^) + a^ea^] (/3a + a^0a) = . 
Regardless of which factor vanishes first, the solution to this last differential equation takes the form 



where C and x ^-re two integration constants and if the first factor of the differential equation vanishes then 

/3fc + 4<6'CX = . 

With A = (3a^) / (2/3), the equation corresponding to Tig yields 

(a A - APA") [a^ {aA^ ~~ 4f3A'^) + 4/3 (/c + a^)] = 

and is therefore autom atica lly satisfied. 

If the first factor in (3.24) vanishes then we have 



p= -p = Q 



The Hubble parameter for this case is 



H 



2 k aS,0 



On the other hand if the second factor in (3.24) vanishes then we have the equation of state 



1 A'(Vi) 1 



[A{Q)r 



where A (y) is given by (3.21). The corresponding expression for the Hubble parameter is 



H 



with a(t) as given in (3.25). In both cases we have 



d a^9 
a " /T 

We notice that if —a£,9/P is positive then the universe is accelerating otherwise a (t) oscillates in time. 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



IV. DISCUSSION 



Among the solutions found in t his an alyses, special attention should be paid to that described by a time dependent 
energy density and pressure as in ( 3.30 ). However, the equations of states p — — p/3 cannot correspond to an ordinary 
dust as the pressure p is negative. Let us suppose that the only matter present on the brane is a self interacting scalar 
field with energy-momentum tensor 



5.0a,</)-ig„ (9'= (/.9fc(/. + V^ ((/))) 



(4.1) 



where gij is a Friedmann-Robertson- Walker metric on the brane. If the field (f) depends only on time then one has 



"-t n — X 



Hy) 



P = Tl = r| = T| = 1 



- V 



Hv) 



(4.2) 
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The field (p is subject to the equations of motion 



with V being the derivative of the potential V with respect to 0. 
The equation of state p ~ —p/2> leads to 



(4.3) 



V = A4> 



(4.4) 



while the equation of motion of the scalar field yields 



Combining these last two equations results in 



3^0 = -V 
a 2 



(4.5) 



+ - = 

6 a 



(4.6) 



Therefore (j) = v/a, for some integration constant v. Using the expression of a (t) as given in (3.25), we find that 



<{t) 



: arctan 



■ exp 



(4.7) 



where w is an integration constant. In order to be able to express the potential y as a fonction of cj), it is convenient 
to rewrite this last equation in the form 



exp 




tan 



1 



TJj) 



(4.8) 



Since the potential is given hyV ^ = 2w2/a2, we find that 



y(0) = — sin^ 
2CX 



/3 



(4.9) 



This potential has an infinite number of minima. 

As it is known, cosmological scenarios with negative pressure are used in explaining the current accelaration of our 
universe (see for example [ |5l[ | ) . The scalar field providing this pressure is known as quintessence [ p2| . Similar models 
are also constructed in the context of brane world cosmology [^-p6| . It seems that wc have found here another model 
for quintessence where the scalar field is governed by a very simple potential. 
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